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Abstract: If new particles are discovered at the LHC, it will be important to determine 
their spins in as model-independent a way as possible. We consider the case, commonly 
encountered in models of physics beyond the Standard Model, of a new scalar or fermion 
D decaying sequentially into other new particles C,B,A via the decay chain D — > Cq, 
C — ► Bl near , B — ► Al iar , Z near and / far being opposite-sign same-flavour charged leptons and 
A being invisible. We compute the observable 2- and 3-particle invariant mass distributions 
for all possible spin assignments of the new particles, and discuss their distinguishability 
using a quantitative measure known as the Kullback-Leibler distance. 
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1. Introduction 

The discovery of new physics at the TeV scale will be a principal objective of experiments 
at the soon-to-be completed Large Hadron Collider (LHC). In most scenarios for physics 
beyond the Standard Model (BSM), new strongly- interacting particles will be observed 
if the collision energy and luminosity are sufficiently high. Typically these particles are 
expected to decay weakly into cascades of Standard Model particles and, possibly, a stable 
or metastable lightest new particle. In supersymmetry (SUSY) with R-parity, for example, 
produced squarks will decay into quarks and, depending on the mass spectrum, leptons 
and/or vector or Higgs bosons and the lightest supersymmetric particle, most often an 
unobservable neutralino. 

Until recently, discussions of new physics at the LHC tended to concentrate on how to 
determine the free parameters of a particular, usually supersymmetric, model. Now that 
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the completion of the LHC is near, however, there is increasing interest in more model- 
independent approaches. Part of this interest arises from the realisation that there are 
BSM scenarios in which the spins of produced particles differ from those expected on the 
basis of supersymmetry. Therefore one needs to consider more generally the ways in which 
the spins of new particles can be determined from their decay chains. 

In the present paper we assume that a particular chain of decays that is common in 
various models, starting from a new scalar or fermionic quark, has been identified and that 
all the masses of the new particles involved in it are known. We then study the extent to 
which decay correlations, manifest in the invariant mass distributions of combinations of 
observable decay products, would enable one to distinguish between the different possible 
spin assignments of the new particles. This paper is thus an extension of earlier work in 
which the SUSY decay correlations were compared with uncorrelated phase space [1, 2] or 
with those of a model that has universal extra dimensions (UED) 1 [4, 5, 6, 7, 8, 9]. 

In the following section we present the decay chains to be considered, and the possible 
assignments of new particle spins and observed particle chiralities. In section 3 we present 
our new results on the corresponding decay correlations. We derive simple analytical 
formulae for the correlation coefficients in terms of the masses in the decay chain, which 
should be of general use whatever the mass spectrum might turn out to be. These extend 
those already given in [4, 10]. We show graphical results for two specific mass scenarios, 
one considered more probable in SUSY and the other in a UED model. In both mass 
scenarios we compare the correlations predicted by all the possible spin assignments. 

As was emphasised by Barr [1], the observability of some of the correlations depends 
on the fact that the LHC is a proton-proton collider, so that scalar or fermionic quarks are 
produced somewhat more copiously than their antiparticles. For the purposes of illustra- 
tion, we take the quark fraction to be f q = 0.7, as was the case for the SUSY and UED 
models studied in ref. [4]. 

In section 4 we consider the extent to which the different spin assignments can be 
distinguished on the basis of a quantitative measure known as the Kullback-Leibler dis- 
tance [11]. This allows us to establish a lower limit on the number of events needed to 
discriminate between any two spin assignments at a given level of confidence. Experi- 
mental effects such as resolution and backgrounds will of course increase the number of 
events needed. However, these are dependent on details of the detector and analysis, and 
so we do not consider them here. The lower limits we compute, being independent of such 
details, can be used to assess whether or not discrimination between two particular spin 
assignments is possible even in principle with a given amount of data. 

We perform analyses of this type on all the observable invariant mass distributions 
separately, and also a combined analysis of the full three-dimensional phase space distribu- 
tion. Our results and conclusions are summarized in section 5. The more lengthy formulae 
are consigned to the appendices. 



lr The similarities between SUSY and UED at hadron colliders were first pointed out in ref. [3]. 
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2. The decay chain 
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Figure 1: The decay chain under consideration. 

We will consider the cascade decay of a heavy colour-triplet scalar or fermion D, of the 
form D — ► Cq, C — > £Z near , 5 -► ,4/ far (figure 1). The decay products are fixed as being a 
quark jet, a pair of opposite-sign-same-flavour (OSSF) leptons, and a stable or long-lived 
massive new particle A. We assume that the masses of the four unknown particles A, B, C 
and D have already been measured (see [12, 13] for example, where edge analysis is used). 
All possible spin configurations in the decay chain are listed in table 1. 
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Table 1: Possible spin configurations in the decay chain (figure 1). 

These 6 chains will be labelled SFSF, FVFV, FSFS, FVFS, FSFV and SFVF respec- 
tively in what follows. Note that SFSF and FVFV are the SUSY and UED cases. 

For fixed spin assignment, there are two possible angular distributions within the chain 
as the quark and near lepton can have either the same or opposite helicity. We will follow 
the conventions of [1] and label these 

• Process 1: {q, / ncar , Z far } = {q L M,lt) or fe>*L>*Z} or UlJrJr} or {q L J-,l+}; 

. Process 2: {q, / ncar , Z far } = {<Zl,£,^} or or {q L ,l^} or {q L ,l+,l R }. 

Note that in some of the processes above (FSFS and FSFV), spin information between the 
quark and near lepton is lost as they are joined by a scalar. For these chains, processes 1 
and 2 give the same distributions. 

Treating the propagators of the unstable particles in the zero-width approximation 
and neglecting all Standard Model particle masses, we can express the matrix elements 
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for these processes in terms of the mass of A, B, C, D and the three two-particle invariant 
masses of the quark plus near lepton, the quark plus far lepton, and the dilepton. It will 
be convenient, as in [4], to introduce the mass ratios 

x = m 2 c /m 2 D , y = m 2 B /m 2 c , z = m 2 A /m 2 B , (2.1) 

so that < x,y,z < 1. The resulting formulae for the full spin correlations are given in 
appendix A. 

To distinguish between the spin assignments, we integrate over two of the indepen- 
dent variables and compare the predictions for the observable invariant mass distributions. 
Throughout, we will show graphical results for two mass spectra of the unknown particles 
A,B,C,D. The first (I) is the MSSM Snowmass point SPS la [14], with fairly widely- 
spaced masses typical of SUSY. The relevant particles and their masses at this point are 
listed in table 2. 
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Table 2: Mass Spectrum I in GeV: Snowmass point SPS la. 

The second mass spectrum (II), with more nearly degenerate masses considered more 
likely in a UED type scenario, is given in table 3, where now the particles involved are 
Kaluza-Klein excitations of Standard Model particles. This UED spectrum was calculated 
using the formulae for the radiative corrections given in [15] with R- 1 = 800 GeV and 
AR = 20. Notice that in this scenario particle C decays into left-handed leptons, whereas 
spectrum I involves right-handed leptons, as is the case for MSSM point SPS la. 
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Table 3: Mass Spectrum II in GeV: Calculated in UED with R- 1 = 800 GeV. 

3. Invariant mass distributions 
3.1 Dilepton mass distributions 

The dilepton mass, run, is the same in processes 1 and 2 and is also relatively easy to 
measure, making it a potentially powerful tool. It depends only on the B decay angle, 
defined as the angle 9 between the two leptons in the B rest frame, through: 

m 2 u = \x{\ - y)(l - z)(l - cos0)m 2 D . (3.1) 
We define therefore the rescaled dilepton invariant mass 

fha = m n l (mn)max = sin(6>/2). (3.2) 
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Figure 2 shows the dilepton mass distribution, dP/dm^, as a function of m 2 > for the 6 
decay chains under consideration for mass spectra I and II. The analytical equations for the 
functions are in appendix B.l. Figures 2-12 are plotted as functions of m 2 , as opposed to 
functions of fh as was done in [4] , as this makes it easier to see the functional dependence. 
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(a) (b) 
Figure 2: Dilepton mass distributions for (a) mass spectrum I and (b) mass spectrum II. 

We see from figure 2 that, as found in [4], the SFSF (SUSY) and FVFV (UED) decay 
chains would be very hard to distinguish on the basis of the dilepton distribution. On the 
other hand the FSFS and FVFS cases, where one or both of the UED vector particles is 
replaced by a scalar, are characteristically different, as is the chain in which one SUSY 
scalar is replaced by a vector (SFVF). Here and in the discussion of subsequent plots, we 
shall quantify these initial qualitative observations in section 4. 

3.2 Quark and near lepton mass distributions 

The quark and near lepton distribution is not experimentally observable as the near and 
far leptons cannot be distinguished. We can however measure jet mass distributions, 
as pointed out in [1]. In order to compute these, we must first calculate the near and far 
distributions. These are then combined in section 3.4. 

The quark and near lepton invariant mass, rn^ r , is given in terms of the angle between 
the two particles, 9*, in the rest frame of particle C: 

KT) 2 = \il - x)(l - y)(l - co S 6*)m 2 D . (3.3) 
We then define the rescaled invariant mass 

m°r = <r7«r)max = ^72). (3.4) 

Figure 3 shows the quark and near lepton mass distribution, dP/d(m^ ear ) 2 , in process 1 
as a function of (wij\ ear ) 2 for mass spectra I and II. Figure 4 shows the same thing for process 
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Figure 3: Quark and near lepton mass distributions for process 1 with (a) mass spectrum I and 
(b) mass spectrum II. 
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Figure 4: Quark and near lepton mass distributions for process 2 with (a) mass spectrum I and 
(b) mass spectrum II. 



2. These are reflections of the distributions for process 1 about the point (fn^f ar ) 2 = 1/2. 
The analytical equations for the functions are in appendix B.2. 

Note that the distributions for FVFV and FVFS and for FSFS and FSFV are the same 
as the differences in the chains only occur at the last vertex, which has no effect on the 
quark and near lepton distribution. Furthermore the distribution for FSFS and FSFV is 
flat since the quark and near lepton are connected by a scalar in these cases. The prospects 
for distinguishing the remaining cases look quite favourable. However, as remarked above, 
we must first take account of the contribution of the far lepton. 
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3.3 Quark and far lepton mass distributions 

The quark and far lepton mass is a more complicated expression than the other two. It is 
a function of all the decay angles in the chain: 

(mlrf = ~{l-x)(l-z) [(l + y)(l-cos0*cos0) 

+ (1 — y)(cos6* — cos#) — 2- v /ysin6>* sin#cos</>] m 2 D , (3-5) 



where 9* and 9 are as before and <fi is the angle between the ql nea,r and dilepton planes, in 
the rest frame of B. The maximum value is at 9* = and 9 = it giving 



Kt r )max = {1 ~ X)(l - Z )m 2 D , 



(3.6) 



so we define the rescaled mass to be 

1 



<f = KfVKf) 



fan 



-.far' 



max " g [(l + y)(l -COS0*COS0) + 



1 -y) (cos 9* - cos 9) -2y/y sin 0* sin 9 cos 0] 5 . (3.7) 

far 



Figure 5 shows the quark and far lepton mass distribution, dP/d(m^ T ) 2 , in process 1 



as a function of (m^) 2 for mass spectra I and II. Figure 6 shows the same thing for process 
2. The analytical equations for the functions are in appendix B.3. 
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Figure 5: Quark and far lepton mass distributions for process 1 with (a) mass spectrum I and (b) 
mass spectrum II. 



It can be seen from figures 5 and 6, and from the formulae in appendix B.3, that in 
many cases the dependence on the quark plus far lepton mass is absent or weak throughout 
the region < (mjf) 3 < y. The main exception is SFVF, where the far lepton is connected 
to the chain via a vector particle. At higher values of fr&, the discrimination is marginally 
better, which is naturally more useful for spectra with smaller values of y (mass spectrum 
I). It is clear, however, that the contribution from the far lepton will degrade the power of 
the quark plus lepton distribution to distinguish between spin assignments. 
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Figure 6: Quark and far lepton mass distributions for process 2 with (a) mass spectrum I and (b) 
mass spectrum II. 

3.4 Observable quark-lepton mass distribution 

We now consider the jet + lepton combinations, jl^, as first used in [1]. As in [4], we 
assume that the jet and lepton are indeed decay products from process 1 or process 2. 
Then the jl mass distribution for a given lepton charge receives near- and far-lepton 
contributions from both the corresponding process and the charge conjugate of the other 
process. In other words, we have 



dP 1 



dm,ji+ 2 



/ dP 2 , dPi_\ / dPi , dP 2 

JQ J near , far I ' J Q \ j near ' 



dm» ear dm^J Jq \dm*f ai dmjf 



(3.8) 



where f q and fq are the quark and antiquark fractions in the selected event sample. The 
factor of one-half enters because Pi 2 are both normalized to unity. Similarly 



dP 1 



dniji- 2 



f ( dPl dP2 1 + f f dP2 dPl 

/</^ dm ncar dm^ ) /9 "^dm^ ear dmj^ 



(3.9) 



This has assumed that both of the leptons are left-handed - in the case of right-handed 
leptons, the expressions for jl + and jl" are interchanged. 

In [4], the Herwig event generator [16, 17] was used to numerically calculate f q in 
the cases of SUSY and UED for two different mass spectra. Despite the differences in the 
models, both gave f q ~ 0.7 for both mass spectra. This is therefore the value that we will 
take for all of our models. 

Figures 7 and 8 show the resulting jet+lepton mass distributions for mass spectra I 
and II, where we have normalised fhji to the maximum observable mass in each case. 

The resulting charge asymmetry, 

_ dP/dm jl+ - dP/dm j7 _ 

dP/dm j i+ + dP/dm j i^ 1 j 
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Figure 7: jet+/ + mass distribution for (a) mass spectrum I and (b) mass spectrum II. 
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Figure 8: jet+Z mass distribution for (a) mass spectrum I and (b) mass spectrum II. 



is plotted in figure 9. There is a reversal of sign between the predictions for the two mass 
spectra, since, as discussed in Section 2, the SUSY mass (I) plots assume right-handed 
leptons, while the UED mass (II) plots assume left. 

It is clear from figures 7 and 8 that high statistics would be required to distinguish 
spin assignments on the basis of the jet plus lepton distributions. This will be seen more 
quantitatively in section 4. On the other hand, the asymmetries in figure 9 show strik- 
ing characteristic differences that would provide convincing confirmation of a model once 
sufficient events had been accumulated. 
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Figure 9: Asymmetry A for (a) mass spectrum I and (b) mass spectrum II. 



3.5 Quark, near and far lepton mass distribution 

The unsealed ql nea,r l iar mass-squared is the sum of the three pairwise unsealed masses- 
squared already considered: 



m 



KT) 2 



iff x < yz, 
iff z < xy, 
iff y < xz, 
otherwise. 



(3.12) 



d - y)(l - z){m u Y + (1 - x)(l - y){m^f + (1 - x)(l - z)(m%) 2 . 

(3.11) 

From eqs. (3.2), (3.4) and (3.7), m 2 u is a complicated function of the three angles 9,9* 
and 4>. The maximum value of this depends upon the relative values of x,y and z. From 
[12, 13], 

m D( 1 ~ x )( l ~ V z ) 

m 2 D (l - xy)(l - z) 
m 2 D (l - xz)(l - y) 
m* D (l - ^/xyl) 2 

Both mass spectra I and II have x < yz; however, in order to remain independent of the 
relative values of x, y and z, we shall no longer define fh q u to lie between and 1. We shall 
instead simply define 

m qll = m qll /m D , (3.13) 

so that it is still dimensionless. 

Figure 10 shows the qll mass distribution, dP/d(m q u) 2 , for process 1 for mass spectra 

1 and II, while Figure 11 shows the same thing for process 2. The analytical equations are 
discussed in appendix B.4. 

In an experimental situation we will be unable to distinguish between processes 1 and 

2 and will instead see a jet+dilepton distribution that is an equal mixture of the two. The 
asymmetry between particle and antiparticle production at a pp collider does not help here, 
because the decays of the unstable particles in the chain into charge conjugate modes must 
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(a) (b) 

Figure 10: qll mass distributions for process 1 with (a) mass spectrum I and (b) mass spectrum 
II. 




(a) (b) 

Figure 11: qll mass distributions for process 2 with (a) mass spectrum I and (b) mass spectrum 
2. 

be equal. Figure 12 shows the combined jll plots for mass spectra I and II. We see that 
there is hope of distinguishing between some cases in the intermediate region of invariant 
mass. 

4. Model discrimination 

It is natural to ask whether from real experimental data it would be possible to determine 
the underlying set of particle spins present in Nature, given only access to one or more of 
the invariant mass distributions of section 3. In particular it would be useful to know how 
many events an experiment would need to identify from one of these decay-chains, in order 
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Figure 12: Jct+dilcpton mass distributions for Processes 1 and 2 combined for (a) mass spectrum 
I and (b) mass spectrum II. 

to favour the assignment of any one spin configuration over any other at a fixed level of 
confidence. We calculate such a number for each pair of spin configurations and for each 
distribution (in isolation from the others) assuming a "perfect" detector with infinite ac- 
ceptance. As real detectors are imperfect and will tend to smear the distributions (thereby 
losing information) the number of events we calculate should be seen as representing a 
lower bound on the number of events that a real detector would need to identify. 

To calculate the number of events ./V needed to disfavour a spin configuration S relative 
to an alternative spin configuration T (which is assumed to be the actual one generating 
the observed events) we solve the following equation for TV 

1 p(S\N events from T) 

— = — — ! -, (4-1) 

R p(T\N events from T) ' 

in which R is the factor by which configuration S is to be disfavoured with respect to T. 
In this note, R was taken to be 1000. Note that the numerator and the denominator of 
the right-hand-side of eq. (4.1) do not transform into each other under the interchange of 
S and T as a result of T's special role as the true configuration chosen in Nature. This 
asymmetry is deliberate, reflecting the fact that there will only ever be one underlying 
distribution which generates the observed events, regardless of the questions which may 
then be asked of the nature of those events. 2 

4.1 One-dimensional analysis 

If we characterise the U N events from T" by the ./V values of a particular invariant mass, 

(T) 

Wb i (for i £ 1,2, N), that are observed in those events, then by using Bayes' Theorem 

2 It would be possible to write p(S\N events from U)/p(T\N events from U) in place of the right-hand 
side of eq. (4.1) by introducing a third model U representing the actual production process. However, the 
large number of choices which could be made for U (beyond the already considered choice of T) do not 
suggest that this course of action would be appropriate in the context of this paper. 
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we can rewrite eq. (4.1) as 

1 p(S) p(N events from T\S) 
R ~ p(T) p(N events from T\T) 

P{S) T^ =l p{m { P\S) 
P^Ul lP (mP\T) 

P(T) \ti \p( m P\T)J) 

In the limit N 3> 1 the sum in eq. (4.2) may be approximated by an integral over the 
allowed masses weighted according to how prevalent they are: 

s"$->("jC*(^)**H 

The integral in eq. (4.3) is one that frequently arises in comparisons of distributions [11] 
and is called the Kullback-Leibler distance 3 of S from T, which we shall denote by: 

KL(T,S)= [ log ( P ^ T 2 ) p{m\T)dm. (4.4) 



K p{m\S) 

Drawing all these threads together we may rearrange eq. (4.3) to give: 

logtf + logfg 

N r (4.5) 

KL(T,S) [ ' 

in the limit of large N. Note that the ratio of the prior probabilities for S and T is present 
in eq. (4.5) as expected - strong prior evidence for S over T should lead to an increase in 
the number of events N from one of these chains needed to discredit S. For the numbers 
presented in tables 4-7 and discussed in the following section, however, we assign equal 
prior probabilities to S and T, thereby removing all prior dependence from eq. (4.5). This 
choice, in effect, looks at the tests in isolation from any pre-existing evidence that might 
favour S over T or vice versa. 

Note that the Kullback-Leibler distance is invariant under diffeomorphisms m — ► f(m) 
of the distributed variable. This means in particular that the number of events iV calculated 
in eq. (4.5) does not depend on whether the distributions are considered to be functions of 
masses or of masses-squared - only the intrinsic information content of the distributions is 
measured. 

4.2 Three-dimensional analysis 

To extract the most information from the data we should compare the predictions of 
different spin assigments with the full probability distribution in the three-dimensional 



3 Note that the Kullback-Leibler distance is not symmetric and so does not define a distance in the usual 
sense. It is however always non-negative, and is only equal to zero when the two distributions are identical. 
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space of mu, m ji+ an d rriji-. The ambiguity between near and far leptons means that this 
given by 

P(mu, m jl+ , rriji-) = ^f q [P 2 (m H ,m j7 +,m j7 -) + P 1 (mu,m jl -,m j i+)] 

+ 7;fq[ p i( m iii m ji+i m ji-) + p 2(mii,m j i-,m j i+)] , (4.6) 

where we use Pi^i^iu w-"f ar , rrJ-f) on the right-hand side, assuming both leptons are left- 
handed, otherwise f q and f q are interchanged. 

Instead of trying to evaluate the three-dimensional generalization of the integral in 
eq. (4.4) analytically, it is convenient to perform a Monte Carlo integration. If we generate 
m iu wi°f ar and m f jf according to phase space, the weight to be assigned to the configuration 

pear = jfar = { - j g 

P + -(m lh m«r, mfl) = \ [f q P 2 (m u , m°f ar , mf{) + /^K.m^, mf)] (4.7) 
while that for Z near = l~ , Z far = 1+ is 

P- + (mu,m*r,™% r ) = \ [f q Pi(m u ,m«r,™f) + /^(%, < ar , m^)] . (4.8) 

In the former case, since the distinction between Z ncar and / far is lost in the data (except 
when interchanging them gives a point outside phase space), we must use eq. (4.6) with 
1+ = / ncar , l~ = Z far in the logarithmic factor of the KL-distance, i.e. the contribution is 



' P + -{mii,m*f ar ,m^\T) + P_ + (m u ,m!j? ,m^ ear |T) . , ,, u . 

l0g 1 P+_(m„ , mJT , m*f\S) + P_ + (m H> mf, rojHS) ) P+ ~ ( m »' m /l T ) ■ ( 4 - 9 ) 

Similarly from the configuration Z ncar = Z", Z far = Z + we get the contribution 



/ r P- + (m^mgr,m^|r) + P + -(m M ,m^,mgr|r) > \ _ fa 

log ; c , ■ ; i — i — r P-+\Tnu, m -\ r , m?r \T ) . 

& \P- + {m ll ,m^ 3I ,m^\S) + P + -(m ll ,m^,rn^ ar \S)J +V 31 ]U 

(4.10) 

Denoting the sum of these two contributions at the ith phase space point by KLj (T,S), 
and summing over M such points, we have as M — > oo 

MlogP 



E,KL 4 (T,5) 



N , (4.11) 



which is the Monte Carlo equivalent of eq. (4.5) when the prior probabilities of S and T 
are taken to be equal. Results for R = 1000 and M = 5 x 10 7 are shown in table 8 and 
discussed in the following section. 



5. Discussion and conclusions 

The results of applying the above one-dimensional analysis to the observable dilepton, 
jet+lepton and combined jet+dilepton invariant mass distributions separately are presented 
in tables 4-7, while the results of the three-dimensional analysis are shown in table 8. 
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(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 




SFSF FVFV FSFS FVFS FSFV SFVF 


SFSF 


oo 


60486 


23 


148 


15608 


66 


SFSF 


oo 


3353 


23 


304 427 80 


FVFV 


60622 


oo 


22 


164 


6866 


62 


FVFV 


3361 


oo 


27 


179 232 113 


FSFS 


36 


34 


oo 


16 


39 


266 


FSFS 


36 


44 


oo 


20 22 208 


FVFS 


156 


173 


11 


oo 


130 


24 


FVFS 


313 


184 


14 


oo 13077 35 


FSFV 


15600 


6864 


25 


122 


oo 


76 


FSFV 


436 


236 


15 


12957 oo 39 


SFVF 


78 


73 


187 


27 


90 


oo 


SFVF 


89 


126 


134 


38 42 oo 



Table 4: The number of events needed to disfavour the column model with respect to the row 
model by a factor of 0.001, assuming the data to come from the row model, for the m\ l distribution 
(a) mass spectrum I (figure 2a) and (b) mass spectrum II (figure 2b). 



(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


(b) 


SFSF FVFV FSFS FVFS FSFV SFVF 


SFSF 


oo 


1059 


205 


1524 


758 


727 


SFSF 


oo 


3006 958 


6874 


761 


1280 


FVFV 


1090 


oo 


404 


3256 


4363 


1746 


FVFV 


2961 


oo 4427 


1685 


2749 


3761 


FSFS 


278 


554 


oo 


418 


741 


870 


FSFS 


914 


4201 oo 


743 


9874 


4877 


FVFS 


1605 


3242 


345 


oo 


1256 


2365 


FVFS 


6716 


1699 752 


oo 


656 


1306 


FSFV 


749 


4207 


507 


1212 


oo 


1803 


FSFV 


720 


2666 10279 


649 


oo 


4138 


SFVF 


813 


1821 


751 


2415 


1888 


oo 


SFVF 


1141 


3517 5269 


1276 


4259 


oo 



Table 5: As in table 4, for the fh^ l+ distribution, (a) mass spectrum I (figure 7a) and (b) mass 
spectrum II (figure 7b). 



(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


(b) 


SFSF FVFV FSFS FVFS FSFV SFVF 


SFSF 


oo 


1058 


505 


769 


816 


619 


SFSF 


oo 


3037 689 


8633 


925 


967 


FVFV 


1090 


oo 


541 


5878 


4821 


445 


FVFV 


2985 


oo 2271 


1431 


4368 


2527 


FSFS 


565 


714 


oo 


1032 


741 


2183 


FSFS 


707 


2297 oo 


526 


9874 


5004 


FVFS 


799 


6435 


882 


oo 


2742 


510 


FVFS 


8392 


1450 525 


oo 


653 


843 


FSFV 


806 


4641 


507 


2451 


oo 


413 


FSFV 


924 


4287 10279 


640 


oo 


4036 


SFVF 


692 


541 


2272 


576 


521 


oo 


SFVF 


1047 


2693 5213 


870 


4041 


oo 



Table 6: As in table 4, for the distribution, (a) mass spectrum I (figure 8a) and (b) mass 

spectrum II (figure 8b). 



As expected from figure 2, the dilepton distribution (table 4) distinguishes very poorly 
between the spin assignments SFSF (SUSY) and FVFV (UED), even with perfect data, 
requiring over 60k events for 1:1000 discrimination in the case of mass spectrum I. However, 
this distribution resolves well between FSFS, FVFS, SFVF and other models. 4 

The jet plus lepton distributions (figures 7, 8 and tables 5, 6) are generally more effec- 
tive than the dileptons in distinguishing between SFSF and FVFV, but the discrimination 
is strongly dependent on the mass spectrum. As was emphasised in [4], a 'quasi-degenerate' 
spectrum like II, with large values of the mass ratios x, y and z in eq. (2.1), usually makes 

4 One should bear in mind that when only a few events are involved, the approximation of replacing the 
sum in eq. (4.2) by an integral becomes invalid, and the discrimination may vary widely between particular 
data samples. 
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(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


( b J 


SFSF FVFV FSFS FVFS FSFV SFVF 


SFSF 


oo 


25630 


241 


1040 82589 


476 


SFSF 


oo 


26712 


189 2213 1686 


421 


FVFV 


27315 


oo 


225 


939 14811 


432 


FVFV 


26391 


oo 


224 1323 1073 


545 


FSFS 


224 


204 


oo 


265 252 


2670 


FSFS 


182 


217 


oo 109 101 


1710 


FVFS 


1009 


906 


278 


oo 1095 


504 


FVFS 


2279 


1373 


116 oo 46742 


210 


FSFV 


73158 


13688 


269 


1124 oo 


557 


FSFV 


1749 


1121 


109 47812 oo 


193 


SFVF 


452 


400 


2749 


493 533 


oo 


SFVF 


405 


528 


1712 196 179 


oo 



Table 7: As in table 4, for the distribution for processes 1 and 2 combined, (a) mass spectrum 
I (figure 12a) and (b) mass spectrum II (figure 12b). 



(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


(b) 


SFSF FVFV FSFS FVFS FSFV SFVF 


SFSF 


oo 


455 


21 


47 


348 


55 


SFSF 


oo 


1053 


21 


230 


194 63 


FVFV 


474 


oo 


21 


54 


1387 


55 


FVFV 


1047 


oo 


27 


135 


190 90 


FSFS 


33 


34 


oo 


13 


39 


188 


FSFS 


33 


42 


oo 


19 


22 175 


FVFS 


55 


67 


10 


oo 


54 


19 


FVFS 


242 


140 


13 


oo 


332 33 


FSFV 


341 


1339 


25 


45 


oo 


66 


FSFV 


189 


194 


14 


315 


oo 37 


SFVF 


62 


64 


143 


19 


79 


oo 


SFVF 


66 


95 


118 


35 


41 oo 



Table 8: As in table 4, for the combined three-dimensional distribution, (a) mass spectrum I and 
(b) mass spectrum II. 



it more difficult to resolve between spin assignments. We see however that there are ex- 
ceptions to this trend, e.g. discrimination between FVFV and FVFS or FSFV. 

The jet plus dilepton distribution (figure 12, table 7), like the dileptons alone, proves 
ineffective in resolving between SFSF and FVFV, but discriminates well between FSFS, 
SFVF and most other models. 

The results in table 8 for the three-dimensional analysis show that, as might be ex- 
pected, this method achieves a discrimination that is better than that of a one-dimensional 
analysis applied to any single invariant mass distribution. This could be particularly useful 
in difficult cases like that of distinguishing between SFSF (SUSY) and FVFV (UED). 

We should stress again that the numbers of events in the tables correspond to perfect 
conditions of signal isolation, resolution and detector efficiency. Realistic conditions would 
most likely require much higher numbers. Nevertheless, the results obtained provide a 
guide to the places where attempts to distinguish between spin hypotheses would or would 
not be worthwhile. Independent of the merits of the hypotheses considered here, the 
Kullback-Leibler distance (4.4) is in our view a useful tool for addressing questions of this 
type. 
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A. Spin correlations from matrix elements 



We give here the analytical formulae for the spin correlations in the cascade decays listed 
in table 1. For brevity we write m q f ar as m qn and m^f as m q f here, and omit propagator 
denominators and overall normalization factors. All the distributions presented in this 
paper are normalized to unity and so such factors are irrelevant. SFSF1 denotes scalar- 
fermion-scalar-fermion process 1, etc. The mass ratios x,y,z are defined in eq. (2.1). 

SFSF1: m\ n (A.l) 
SFSF2: (1 - x)(l - y)m 2 D - m 2 qn (A.2) 
FVFV1: x(l - x)y(l - y)(l - z)m% - 2x(l - z)[(y - 2xy - 4z + 2yz)m 2 D 
+4zm 2 qn ]m 2 qn - 2x(l - 2z)[(l - y)m 2 D - 2m 2 qn ]m 2 qf 

+ (1 - 2z)[(l -x){2- y)m 2 D - 2(1 - 2x)m 2 qn ]m 2 (A.3) 

FVFV2: x(l -x)(l- y)(l - z)(y + 8xz)mj) - 8xz(l - z)m qn 

-(1 - 2z)[2x(l - 2x)(l - y)m 2 qf - (1 - x){2 - y)ml]m 2 D 
-2{x{\ - z)[y - 2(1 - 2x){2 - y)z]m 2 D - 2x(l - 2z)m 2 qf 
+ (1 - 2z)m 2 u }m 2 qn (A.4) 
FSFS: x{l - y)(l - z)m 2 D - ml (A.5) 

FVFS1: 2{2xm 2 qf - (1 - 2x)[xy(l - z)m 2 D + m 2 u }}m 2 qn 

+{x(l - y)[(l - x)y(l - z)m 2 D - 2m 2 qf ] + (1 - x){2 - y)m 2 u }m 2 D (A.6) 

FVFS2: -2[xy{\ - z)m 2 D - 2xm 2 qf + m 2 u }m 2 qn - {2x(l - 2x){\ - y)m 2 qf 

-(1 - x)[xy(l - y)(l - z)m 2 D + (2 - y)m 2 ]}m 2 D (A.7) 
FSFV: x{\ - y){\ - z)m 2 D - (1 - 2z)m 2 u (A.8) 

SFVF1: -2x[x{\ - y){2 - z)m 2 D - 2m 2 n ]m 2 qf + 2(1 - x){2x 2 y{\ - y){\ - z)m% 
+ [x(l - 2y)(2 - z)m 2 D - 2m 2 l ]m 2 } + xz[x(l - 2y)(l - z)m 2 D 
-2m 2 ]m 2 qn (A.9) 

SFVF2: -z[x(l - 2y)(l - z)m 2 D - 2m 2 u ]m 2 qn + 2[x(l - y){2 - z)m 2 D - 2m 2 u }m 2 qf 

+ (1 - x)z[x(l - y)(l - z)m 2 D - 2m 2 u ]m 2 D (A. 10) 

B. Analytical formulae for invariant mass distributions 

This section contains all the analytical formulae for the invariant masses for the 6 chains. 
Results for the SFSF and FVFV chains have appeared in the literature before [4, 10], but 
are included here for completeness. Throughout, m 2 represents the relevant rh 2 for that 
subsection. 

B.l Dilepton invariant mass distributions 

The following table (9) contains the dilepton invariant mass distributions, dP/dm 2 ,, for the 
different chains. These are equal in both processes 1 and 2 and have been normalised to 
unit area. 
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Chain 


Processes 1 and 2 


SFSF 


1 


FVFV 


(2+v)(i+2*) (V + 4z + ( 2 - f )(! " 2z ) m2 ) 


FSFS 


2(1 -m 2 ) 


FVFS 


_2_( y + m 2 (2 _ y)) 


FSFV 


T ^(l-m 2 (l-2z)) 


SFVF 


_jy__ {4y + z + 4m 2 (1 _ y{2 -z)-z)- 4m 4 (l - y)(l - z)) 



Table 9: Dilcpton invariant mass distributions. 



B.2 Quark and near lepton invariant mass distributions 

The following tables (10 & 11) contain the analytical forms for the quark and near lepton 
invariant mass distributions, dP/d(m° / ear ) 2 . These have been normalised to unit area. 



Chain 


Process 1 


SFSF 


2m 2 


FVFV 


(i+2x)(2+ W ) Iv + 4 ( x y + xyW 4(1 x)(i y )m 4 ] 


FSFS 


1 


FVFS 


(i+2x)(2+ y ) \y + 4(1 - y + xy)m 2 - 4(1 - x)(l - y)m A ] 


FSFV 


1 


SFVF 


2(2j/+(l-2y) m ^) 
l+2y 



Table 10: Quark and near lepton invariant mass distributions for process 1. 



Chain 


Process 2 


SFSF 


2(1 - m 2 ) 


FVFV 


d+2x)(2+y) I 4x + y + 4 ( x 2x v + xyW 4 (! ^)(! y)™ 4 ] 


FSFS 


i 


FVFS 


(i+2xH2+y) t 4a: + 2/ + 4(1 - 2x - y + xy)m 2 - 4(1 - x)(l - y)m 4 ] 


FSFV 


1 


SFVF 


2(l+(-l+2j/)m a ) 
l+2y 



Table 11: Quark and near lepton invariant mass distributions for process 2. 



It is to be expected that the FVFV and FVFS distributions and the FSFS and FSFV 
distributions match as the chains are identical up to the first two vertices in the chain 
and the quark and near lepton invariant mass is unaffected by what happens at the third 
vertex. 

B.3 Quark and far lepton invariant mass distributions 

The following gives the quark and far lepton invariant mass distributions for the different 
processes, dP/d(m^ r ) 2 , normalised to unit area. The index 1 or 2 indicates Process 1 or 2 
respectively. 
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SFSF 



dPi _ -2 

dm, 2 (1 — y) 2 



dPo 



FVFV 



dm 2 (1 — y) 2 



6 



(1 — y + logy) < to 2 < y 

(1 — m 2 + logm 2 ) y < m 2 < 1 

(1-y + ylogy) < m 2 < y 

(1 — m 2 + y log to 2 ) y < m 2 < 1 



dPi = 

dm 2 (1 + 2x)(2 + y)(l + 2z)(l - y) 2 " 

' (1 - y)[4x - y + 2(2 + 3y - 2x(5 + y))z 

-4m 2 (2 - 3x)(l - 2z)\ - [y(l - 2(4 + + Ax{2z - y(l - 4z)) 
+4m 2 (l + y - x(2 + y))(l - 2z)\ logy 



(B.l) 



(B.2) 



< m 2 < y 
(B.3) 



(1 - m 2 )[4a;(l + 2y - 5z - <oyz) - 5y + 2(2 + 9y)z 

-4m 2 (l - x)(l - z)\ - [y{\ - 2(4 + y)z) + 4x(2z - y(l - Az)) 

+4m 2 (l + y - x(2 + y))(l - 2z)] logm 2 y < m? < 1 



dP 2 _ 6 

dm 2 ~ (l + 2x)(2 + y)(l + 2z)(l-y) 2 X 

' (1 - y)[-y + 2(2 + 2x(l - y) + 3y)z - 4m 2 (2 - x)(l - 2z)\ 

-[y(l - 2(4 + y)z) + Am 2 (l + (1 - x)y)(l - 2z)} logy < m 2 < y 

(B.4) 

(1 - m 2 )[4(l + x)z - y{5 - 18z + 8xz) - 4m 2 (l - x)(l - z)\ 

-[y(l - 2(4 + y)z) +4m 2 (l + (1 - x)y)(l - 2z)]logm 2 y < m 2 < 1 



FSFS 



dPi,2 = -2 
dm 2 (1 — y) 2 



(1 - y + log y) < m 2 < y 
(1 — m 2 + log to 2 ) y < to 2 < 1 



(B.5) 



FVFS 



dPi _ 6 

dro 2 ~ (l + 2x)(2 + y)(l-y) s 



(1 -y)(4x-y- 4m 2 (2 - 3x)) 
+ [(-l + 4x)y 

+4to 2 (2x — 1 - (1 - x)y)] log y < to 2 < y 

(B.6) 

(1 — to 2 )(4x — 5y + 8xy — 4m 2 (l — x)) 
+ [(-1 + 4x)y 

+4to 2 (2x — 1 — (1 — x)y)] log to 2 y < to, 2 < 1 
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dP 2 6 

dm 2 ~ {l + 2x){2 + y){l-yf 

FSFV 

dPi, 2 -2 



{ {l-y){-y-Am 2 {2-x)) 

— (y + 4m 2 (l + (1 — logy < m 2 < y 



(B.7) 



(1 - m 2 )(-5y - 4m, 2 (l - x)) 

— (y + 4m 2 (l + (1 — x)y)) logm 2 y < m 2 < 1 



dm 2 (l + 2z)(l-y) 2 
SFVF 

6 



(l-y)(l-22) + (l-2y2)logi/ < m 2 < y 



(B. 



(1 - m 2 )(l - 2z) + (1 - 2yz)\ogm 2 y < vn? < 1 



dPi 



dm 2 (l + 2y)(2 + z)(l-y) 2 



' (1 - y){2 - 3z - 2y(l + z) + 4m 2 (l - 2z)) 

— (2 + Ayz — 4m 2 (l — z — yz))logy < m 2 < y 

(1 - m 2 )(2 -3z- 8yz + 2m 2 (l - z)) 

— (z + 4yz — 4m 2 (l — z — yz)) logm 2 y < m 2 < 1 

(B.9) 



dP 2 = 6 

dm 2 ~ {l + 2y){2 + z){l-yf 



' (l-y)(z -4m 2 (l-2z)) 

+(yz — 4m 2 (1 — z — yz)) log y < m 2 < y 

(1 - m 2 ){z{l + Ay) - 4m 2 (l - z)) 

+(yz — 4m 2 (1 — z — yz)) logm 2 y < m 2 < 1 

(B.10) 



B.4 Quark, near and far lepton mass distributions 

Due to the complicated nature of the qll mass, these distributions can be lengthy and 
complicated. Included in this section are the ones of manageable length - FSFS, FSFV and 
SFSF (processes 1 and 2). The others are available from the authors on request. 
These distributions are given in terms of 



M-m = min < 



' i(l-^)(^(l-V^) 2 ~rn^+^(l +v ^)^m^) 2 -(l-x-m^) 
x(l-y)(l-z) 



(B.11) 



This is equal to 1 in the region between (1 — y)(l — xz) and (1 — z)(l — xy), which leads 
to the flat section of the distributions observed. These distributions require the positive 
constants Mi to be set such that the distribution integrates to 1 in each case. 

FSFS 



dp 



dm 2 , (l-y)2(l_ 2 )2 



(B.12) 
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(1 - yz) - ^((1 + VP) 2 - (1 - y)(l - z)M+)((l - VP) 2 - (l-2/)(l- 4M+) 
-2(y + 2) log 



'(1 + VP) 2 - (1 - y)(l - z)A4 + - J(l - ^~z) 2 - (1 - y)(l - z)M+ 



2v/p 



FSFV 

dP A/2 



dm 2 , (i_ y) 2 (1 _ z) 2 



(l-2z)(l-yz) (B.13) 



-(1 - 2z)yJ((l + VP) 2 - (1 - y)(l - z)M+)((l - VP) 2 " (1 " " *)M] 
-2(y + z-2z(l+yz)) x 



log 



'(1 + jyzf - (1 - y)(l - z)M+ VP) 2 - (1 - y)(l - z)M+ 



SFSF 



A4+(l - y)(l - z)((l + xz)(l + yz) + 2z(l + xy) - rh 2 qll (l + z)) 



^(1 - vp) 2 - aot - y)(i - *) va + VP) 2 - -Hta - " *) 



(1 - yz)(m 2 Jl - z) - (1 + xz)(l - yz)) 



^(1 - VP) 2 - - m - z)yj{\ + VP) 2 - M+(l - y)(l - z) 

\l + VP) 2 - (1 - y)(l - z)M+ - - VP) 2 - (1 - m - z)M+ 



-4(1 + xy) log 



2VP 



dP 2 dPi 



d ™ 2 « d ™ 2 ;; 



(B.15) 




8AA3 , l\l{l + Jyz) 2 -{l-y){l-z)Mt l - x /(l-^yz) 2 -(l-y)(l-z)Mt l 



(l-y)(l- Z ) 2 VP 



An explicitly normalised expression for an equal mixture of process 1 and process 2 
for the SFSF (supersymmetric) spin configuration may be found in [18]. 
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